ON THE EULERIAN POLYNOMIALS OF TYPE D 

CHAK-ON CHOW 



Abstract. We define and study sub-Eulerian polynomials of type D, which count the el- 
ements of Dn with respect to the number of descents in a refined sense. The recurrence 
relations and exponential generating functions of the sub-Eulerian polynomials are deter- 
mined, by which the solution to a problem of Brenti, concerning the recurrence relation for 
the Eulerian polynomials of type £), is also obtained. 



1. Introduction 

One of the classical polynomials of combinatorial significance is the Eulerian polynomial, 
which enumerates elements of the symmetric group with respect to the number of descents, 
and whose properties are well studied (see, e.g., 0). 

The Eulerian polynomials and their g-generalizations have also been defined for other 
Coxeter families (see, e.g., for g- Eulerian polynomials which interpolate between the type 
A and -B, and between type A and cases). For instance, the basic properties of the 
Eulerian polynomials of type i?, analogous to their type A counterparts, are known. 

On the other hand, the type D theory is not as well developed as the type B one does. Only 
the type D generating functions and Worpitzky identity are known; the type D recurrence 
relation is still missing. Finding such a type D recurrence happens to be a problem not as 
simple as in the type A and B cases. The difficulty lies in that, by imitating the derivations 
of the corresponding type A and B recurrences, the number of descents changes in a peculiar 
way so that some refined Eulerian polynomials are needed in order to capture these changes. 
This leads to the introduction of the sub-Eulerian polynomials, which is the idea central to 
the present work. 

The organization of this paper is as follows. In the next section, we collect the notations 
that will be used in the rest of this work. In section 3, we introduce the sub-Eulerian 
polynomials, which enumerate the elements of the Coxeter group of type Z) in a refined 
sense. In particular, we obtain a recurrence relation involving the type D Eulerian and 
sub-Eulerian polynomials. In section 4, we compute the exponential generating functions 
of the sub-Eulerian polynomials. The generating functions computed enable us to refine 
the recurrence relation obtained in section 3. In the final section, we determine the partial 
differential equation of minimal order which the exponential generating function of type 
D Eulerian polynomials satisfy, and by which a recurrence involving the type D Eulerian 
polynomials and its derivatives only is derived. 
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2. Notations 

Wc collect some definitions, and notations that will be used in the rest of this paper. Let 
S* be a finite set. Denote by the cardinality of S. Denote by ©„ the symmetric group 
of degree n, Bn the hyperoctahedral group of rank n, and Dn the Coxeter group of type 
D of rank n. Let tt = 7ri7r2 • • • 7r„ e Dn, where 7r(i) — iTi, for i = 1,2,... ,n. We say that 
i E [0, n — 1] is a D-descent of tt if tTj > TTj+i, where ttq = — 7r2. We shall drop the type 
designation of descents in the sequel unless circumstances demand the contrary. Denote by 
Des(7r) the descent set of n, and d{n) = # Des(7r) the number of descents of n. Denote by 
Dn^k the Eulerian number of type D, which is defined as the number of elements of with 
k descents. Define the Eulerian polynomial Dn{t) of type D by 

n 

3. SUB-EULERIAN POLYNOMIALS 

The goal of the present work is to study the recurrence relation satisfied by Dn^k- Towards 
this end, we have the symmetric group and the hyperoctahedral group -B„ as our guiding 
examples. Recall that &n (resp., Bn) can be constructed by inserting n (resp., ±n) to the 
elements of &n-i (resp., Bn-i). By studying how the descent number changes, recurrence 
relations for the Eulerian numbers of type A and B are obtained. 

In the case of Z)„, a similar construction is also possible. Let a = g\<J2 ■ ■ -On-i G ^n-i- 
Denote by a the element a\Gi - ■ ■ Gn-\ of -B„-i. It is clear that the map -B„_i Bn-i, 
(J ^ o", is a bijection sending Dn-i onto Bn-i \ Dn-i- To each element a of Dn-i, we can 
associate to a the pair of elements {o", o"} of Bn-i- This can also be realized as the orbit 
of a under the right action of the subgroup (sq) of Bn-i generated by sq = (1 1) (in cycle 
notation). (It is unfortunate that this right action of (sq) does not endow the collection of 
orbits with a group structure because (sq) is not normal.) With this association in place, 
Dn is constructible from by inserting n to a and —n to a. 

For n ^ 2, denote by ^Dn^k (resp., ^'^Dn,k, ^"^Dn^k, and ^'^"^Dn^k) the collection of elements 
of Dn of k D-descents, with descent at position 1 but not at (resp., with descents at both 
positions and 1, with no descent at positions and 1, and with descent at position but 
not at 1). Denote by ^-D„ fc (resp., "'^D^^^, ^^D„ fc, and °'^^D„ fc) the collection of elements 
of Bn \ Dn of k D-descents, with descent at position 1 but not at (resp., with descents at 
both positions and 1, with no descent at positions and 1, and with descent at position 
but not at 1). 

Let (7 = (7i(72 ■ ■ - (Tn ^ Dn be of k descents. It is of interest to see how the pair {cr, a} 
distributes amongst ^Dn,k, ^Dn,k, etc. 
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Lemma 3.1. We have 



Proof. The map o" ^ o" does not alter (72 • • • so that descents beyond the position 1 are 
unaffected. It remains to figure out how descents at positions and 1 change. But the 
following calculations 



(7i > (72 j > (71 + £72 

^cO,ln (7i + (72<0l CTi>(72 I ;^-f=0,ln 

(7i > (72 I L) > (7i + (72 



^p^2r) (7i + (72>0l C7l<C72 I A p ^2 f) 

(7e-L'„,, ^ (7i<(72 0<(7i+(72i^^ 



^^0,^2r) (7i + (72<0l (7i >(72 I ;V ^ 1 f) 

(7i < (72 L) < (7i + (72 ' 



then show that either the descents at positions and 1 are preserved, or the descent at 
position of (7 is mapped to the descent at position 1 of a, and conversely. In any case, 
the number of descents are preserved, and the leading descent structure are given as in the 
lemma, concluding the proof. □ 



We are now ready to look at the insertion process. Let a = a\a2 ■ ■ ■ (7„-i G *Dn-i,k, where 
*Dn-i,k is any of the subcoUections of elements of D^-i of k descents, defined above. The 
letters ±n can be inserted at any of the n positions, numbered consecutively from to n — 1, 
and the insertion proceeds as follows. 

At position 0: 

if (7 e ^Dn-i,k (o- e °'^^Dn-i,k), then 

(1) n(7i(72 • • • (Tn-l e ^-Dn,fc+1, riCTi(72 • • • (7„_i G °'^^-Dn,fc; 

if (7 G O'^L'n-i.ifc {a G "'^L'n-i.fc), then 

(2) n(7i(72 • • • an-1 G ^-D„,fc, nai(72 ■ ■ ■ (7„_i G ^'^'^Dn,k; 
if (7 G ^^A^-Lfe (-7 G ^2Z:'„-i,fc), then 

(3) naia2 ■ ■ ■ (Ja-i G ^ A,fc+i, ^ai(72 • • • a^-i G °'^^-D„,fe+i; 
if (7 G ^'^^Dn-i,k (o- G ^i^n-i.fe), then 

(4) 71(71(72 • • • (7„_i G ^-D„,fc, 72,0-1(72 • • • (7„_i G °'^^£)„,fc+i. 

At position 1: 

if (7 G ^£>„-i,fe {a G ^'^^D^_i^k), then 

(5) (7i7l(72 • • • (7„_i G ^'^Dn,k, 0x7102 ■ ■ ■ (7n-l G ^'^Dn,k+l] 



4 CHAK-ON CHOW 

if fT G °'iZ^„_i,fc (a G °'iZ^„_i,fc), then 

(6) o-ino-2 • • ■ o"n-i e ''^Dn,k-1, ^in(T2 ■ ■ ■ 0-n-i ^ °'^I^„,fc; 
if a G ^"^Dn-i^k e ^^/^n-i.fc), then 

(7) <Tin(j2 ■ • ■ o-„_i G ^^/^n,fc+i, o-incxa ■ ■ ■ cr„_i G °'^£'n,fc+2; 



if (T G "'^^/^n-i.fc ((7 G ^/^n-i.fc), then 



(Tincr2 ■ ■ ■ cr„_i G =*^-D 



At position i [2 ^ i ^ n — 2): 
if cr G *Dn-i ki then 



(9) 



(Ti(T2 ■ ■ ■ aiUai+i ■ ■ ■ Gn-l G < * 



0-1^0-2 ■ • • cr„_i G °'"^-D„,fc+l- 



"L'n^fc if i G Des(a), 
'L'^^fc+i if^^Des(a); 



if a G *Dn-i,k, then 



(10) (Ticr2 ■ ■ ■ (Tiuai+i ■ ■ ■ an-i G 



*Dn,k if i e Des(a; 
'-D„,fc+i if i ^ Des{a] 



At position n — 1: 

(11) 

(12) 



cricr2 • • • CTn-in G *Dn,k', 



^n.k 



Define now the suh-Eulerian numbers Dl^j., -D°'^^, -D^^^, -D^'J"^ by 

^ Dn,kj 

■,ki 



^^2 



^n,k 



The sub-Eulerian numbers defined are obviously refinements of the Eulerian numbers Dn,k-, 
and they offer a more accurate description of the descent distribution of elements of Dn- 



Since the map cr ^ o" is an injection, Lemma |3.1| then says that 



iP'Dn,k — 


D 


J±0,l f) _ 
TT ^n.k 


D 


ih^'Dn.k = 


D 




D 



0,^2 
n,k 1 

0,1 
n,ki 

^2 
n,ki 

1 

n,k' 



Note that, in the insertion process above, the descent numbers can go up by 2 or go down 
by 1, which never occurs in the case of (3„ or B^. However, the introduction of sub-Eulerian 
numbers helps capture these changes. 
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Proposition 3.2. For n ^ 3, the sub-Eulerian numbers satisfy the following recurrence 
relations 



fii) D, 



0,1 
n.k 
,>2 



{2k 
2{k 



-,0,1 



>2 

n-l,fe-l 



D 



0,^2 
i-l,fc' 



l)^:'-M + (2r^-2/c+l)D°'Vl 



-^n-l,fe-2 ^ -^n-l,fc-l' 
,0,1 , n0,>2 



iii) Djt = {2k + 1)D^\, + 2{n-k- l)D^\,_, + Dl_,^, + Z^l^'^^^,^, + D^^,^,,, 



fiv) D 



0,^2 
n.k 



+ (2n - 2A; - l)D 



0,^2 
n-l,A:-l 



^n-1 



,fc-l + -^n-l 



^2 

n-l,fc-l- 



Proof. We only prove (i); the remaining assertions follow from similar considerations. Ele- 
ments of ^Dn^k can be obtained by inserting n 



io a e ^Dn-i,k^i 
o G °''/^n-i,fc 

^ c 0,^2 n 

a e Un-l,k 

a e ^Dn-i,k 
a e ^Dn-i,k 

and by inserting —n 

to cr e ^-D„_i,fc 
a G ^-Dn-i,fe-i 
a G ^-D„_i,fc_i 



as in (|T]) with multiplicity -D^_i 



0,1 

n—l,k 
^2 

n— l./c— 1 
0,>2' 



as m 



m 



D 
D 

^n-l.k 

{k - l)Di_,^, 
[{n -2)-{k- 

^n-l.k 



with multiplicity {k — l)Dl^_^ ^, 

{(n -2)-{k- . 1, 



^n-l,k-l- 



Summing the multiplicities, we finally have 

'n,k 



Dl , = {2k- l)Dl ., + 2{n - k)Dl^,,_, + 



n— l.fc 



D 



^2 

n— l.A;— 1 



D 



0,^2 
n— l,fc' 



as desired. 



□ 



It is convenient to record the sub-Eulerian numbers by generating functions. For n ^ 2, 
define the sub-Eulerian polynomials D\{t), D^^{t), D^'^{t), and D°'^^(t) by 



Di{t) 

D'r.'\t) 

D'^\t) 



k>0 



fc^O 



fc>0 



k>0 



The first few values of D\{t), etc., are given in Tables |T|-§. 

It is clear that the Eulerian polynomial Dn{t), n ^ 2, satisfies 
(13) D^{t) = Dl{t) + D'/{t) + (t) + 
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The Eulerian polynomials of type A and B satisfy certain difference-differential equations. 
The same is true of the sub-Eulerian polynomials. 

Proposition 3.3. For n ^ 3, the sub-Eulerian polynomials satisfy the following difference- 
differential equations: 

(i) Diit) = [2{n - l)t - l]Di_,{t) + 2t{l - t){Dl_,y{t) + D'l,{t) + tDl\{t) + D^l{t), 

(ii) D'/{t) = [{2n-l)t-2]D'l,{t) + 2til-t){D'lJ{t)+tDl^_,{t) + t'Df_,^ 

(iii) (t) = [2{n-2)t + l]Df^,{t) + 2til~t)iDf_,)\t) + Dl^_,it)+t-'D'/^^^^^ 

(iv) = i2n - 3)tD'^^^{t) + 2t{l - t){D',:^l)'{t) + tD^^^it) + D^^,{t) + tDf_,{t). 

Proof. We only prove (i); the remaining assertions follow from similar reasoning. Multiplying 
Proposition p.2| (i) by t^ and summing over k yields 

ran n n 

k=l k=l k=l k=l 

k=l k=l 
= I + II + III + IV + V. 

The left hand side of (p!4D is equal to -D^(t) because D\ q = 0, while terms on the right hand 
side are equal respectively to 

n n 

I = 2tY,kDl_,^,t'-' -Y,Dl_,,kt' = 2t{Dl_,y{t) - Dl_,{t), 

k=l k=l 

n— 1 n—l n— 1 

11= Y.^{n-k- l)D\_,^,t'^' = 2{n - l)t Dl_,^,t' - 2t' kD\_,^,t'-' 
= 2{!i-l)tDi_,{t)~2t\Di,,y{t), 

n—l 



111= E^nV = ^°-lW' 



fc=0 
n-1 



IV = Y.Df-i,kt'^' = tDr_Aty 

k=0 

v= D'^^^ity 

n,0 n— 1,71 n,0 n—l,n n— 1,0 



because D^q = D^.^ = D% = D'^^^\ = D'^^^Iq = D^^'^l = 0. Hence, (i) follows. □ 



In view of the decomposition property (|T3D of the sub-Eulerian polynomials, the following 
is the immediate consequence of Proposition p.3|. 



Corollary 3.4. We have 

Dn{t) = 2ntDn-i{t) + 2t{l - t)D'^^,{t) + (t"! - t)/^°'!i(t) + (1 - tfD^\{t) 

+2{l-t)D'^^^ity 
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Proof. Summing Proposition p.3| (i)-(iv) and using ([131) . □ 

The above recurrence relation will be revisited later when further properties of the sub- 
Eulerian polynomials are found. 
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2 t 

3 3t + 

4 7t + 34t2 + rt^ 

5 15t + 225^2 + 225t3 + 15t^ 

6 31t+ 1196t^ + 3306t=^ + 1196t^ + 31t5 

Table 1. The sub-Eulerian polynomial D\{t) for n = 2, . . . ,6. 



n 




2 




3 


he + e 


4 


ne + 3oe + e 


5 


49^ + 303^ + me + e 


6 


129^2 + 2132t3 + 3030t^ + 468t^ + ^ 



Table 2. The sub-Eulerian polynomial D^^{t) for n = 2, . . . ,6. 



2 1 

3 1 + 5t 

4 1 + 30t + lie 

5 1 + 127t + 303t2 + 49^3 

6 l + 468t + 3030t^ + 2132t^ + 129t^ 

Table 3. The sub-Eulerian polynomial D^'^{t) for n = 2, . . . ,6. 



n 




2 


t 


3 


3t + 3e 


4 


7t + 34^2 + re 


5 


15t + 225^2 + 225^3 + 15t^ 


6 


31t + 1196^2 + 3306t3 + 1196t^ + 3ie 



Table 4. The sub-Eulerian polynomial for n = 2, . . . ,6. 



4. Generating Functions 



Proposition |3.3| offers an efficient way of computing Dl^{t), etc. This section discusses a 
convenient way of recording them, e.g., by their exponential generating functions. Define 
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the exponential generating functions for the sub-Eulerian polynomials by 



^ n 
(15) ^ 

D^\x,t) = J2Dfit)^^ 



n>2 



D'^^\x^t) = Y,D'^\t) 



n>2 



X 

n\ 



Here, we are interested in obtaining closed form expressions for the right hand sides of (|T^. 

Denote by and Bn{t) the n-th Eulerian polynomials of type A and B for 6„ and 

En, respectively. In the case of (3(x,t) = X]n>o ®"(^)'''"/^'' determines fn^k for which 
6„(t)/(l — tY'^^ = J2k>o fnM^ y then multiply by x"'/n\, sum over n, and replace x by 
x(l — t). This same procedure is complicated for D^{x,t), etc., because they are coupled 
via Proposition p.3| (i)-(iv). Also, it is not clear at present what are the right denominators 
Q(t) for the rational generating function D^(t)/Q(t). 

An alternative, which we believe to be new, way of computing them is as follows. Let 
us show how it works in the case of &{x,t). By taking partial derivatives of &{x,t) with 
respect to x and t, and making use of the difference-differential equation which satisfy, 
namely, 

(16) 6„(t) = [{n - l)t + l]6„_i(t) + t{l - t)6Li 

we obtain that & = &{x,t) satisfies the following first order linear partial differential equa- 
tion 

(17) t(t-l)6t + (l-xt)6, = 6, 

which, together with the initial condition (5(0, t) = 1, uniquely determine (5, that is, 

(18) &{x,t)= . 

For completeness, we record the PDE which B{x,t) = J2n^oBn(t)x^/n\ satisfies, namely 

(19) 2t{t-l)Bt+ il-2xt)B^ = {l+t)B. 

The concerned initial condition is -8(0, t) = 1, and the solution to ( |19|) is of course 

(20) Bix,t)='-^-±^. 

We shall not give the details of the proof of (|I^ and (0) but encourage interested readers 
to work them out by imitating the proof of Proposition O below. 

Proposition 4.1. We have 
(i) 2t{t - l)Dl + (1 - 2xt)Dl = -D^ + Z^O'i + W^^ + D^'^^ + xt; 
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(ii) 2t{t - + (1 - 2xt)Dl^^ = tD^ + {t- 2)D°'^ + eo^^ + tD^'^^ ^ ^^2 . 

(iii) 2t{t - l)Df + (1 - 2xt)Df = D^ + t-i/^O'i + (1 - 2t)D^2 ^ ^0,^2 ^ ^ . 

(iv) 2t(t - + (1 - 2xt)Z^p2 _ ^^1 ^ ^0,1 + ^^^2 _ ^^0,^2 ^ 

Proof. We only prove (iii); the remaining assertions follow from similar reasoning. We have 
2t{l-t)Df= 5^2t(l-t)(Df)'(t)^ 

n^2 



n>2 



n^2 ■ n^2 ^ '' 



+ (2t-l)5^D„"^^(t)^-$^DP^(t)^ 
^ — ^ n! — ^ n! 

ra^2 n^2 



n! ^ n\ 

n^2 n^2 

= _ /}2^2^t)a; _ 2a;tDf + (2t - l)/}^^ _ ^0,^2 _ ^1 

from which (iii) follows, because Df^{t) = 1. □ 



The above set of PDEs are readily solved by the method of characteristics . The idea 
is to solve for characteristic x = x{t; A), which satisfies dx/dt = (1 — 2xt)/2t{t — 1) and is 
parametrized by A, where A is the constant of integration; for x and t related by x = x{t\ A), 
the partial derivatives become total derivatives with respect to t, and the PDEs are thus 
reduced to ODEs and the solutions of which solve the given PDEs. (Here, the constants 
of integration of the latter ODEs are arbitrary functions of A, and whose forms can be 
determined by the initial conditions.) 

Theorem 4.2. We have 

f(px(l-t) _ i\2 

(i) D^^^\x,t) = , ^V— = D\x,t), 

V ^ V , ; 2(1 -t)(l -te2^(i-*)) ^ ' 

_ -1 - a;(l - t) + e^(i-*) 



(ii) D^^(x,t) ^ , o n A^ ; 

^ni. ^ t2^(l-*) -t2[l_a;(l_t)]e2-(l-*) 

iii ^ ^ 1 y J ^ 

^ ' ^ ' ' (1 -t)(l -te2^(i-*)) 

Proof. Let x = x{t) be such that dx/dt = (1 — 2xt)/2t{t — 1). The latter linear first order 
ODE is readily solved, e.g., ti/^e^*^^"*-* = A, where A is a constant. On x = x{t), the left 



hand sides of Proposition 4.1 become total derivatives with respect to t. In particular. 
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Proposition [4.1| (iv) is then 

(21) 2t{t - 1)^^ = tD^ + Z^O'i + tD^^ - tZ}0'^2^ 

Differentiating (|21|) with respect to t, using dx/dt = (1 — 2xt)/2t{t — 1), multiplying by 
2t(t — 1), followed by substituting the t-derivatives on the right by the corresponding right 
hand side in Proposition |4.1| , we have that D^'^^ satisfies the following ODE, 

(22) At{t - 1)'^^^^- + - l)(3i - 1)^^^^ - 4^°'^' = 1' 
whose solution (obtained by Mathematica) is 

(23) D°'^'(x,t) 



>2. 1 - t + 4(1 + t)Ci{A) - 8tt'/^C2{A) 



Ait - 1) 

Here, Ci, and C2 are arbitrary functions of A. The initial condition D°'^^(0,t) = implies 
that 

1 - 1 + 4(1 + t)Ci{t'/^) - m^'^C2{t^'^) 

Ait - 1) 

so that 

_ l-t^ + 4(l + t^)Ci(t) 

^^^^^ 8^^ • 

Reporting this back in (123) and replacing A by t^/^e^^^"*) yields that 



_ (e-(^'*) - l)[te-(^-*) + 1 + 4(1 - te-(^-^))Ci(tV^e-(^-^))] 
i24j 1^ 4(t-l)e-(i-*) 

The initial condition D2'^^(0,t) = then yields that 

1 + 4(1 -t)Ci(tV^)) 
A{t - 1) 

from which we deduce that 

t^ + 1 



4(t2 - 1)- 

Substituting Ci{t) back in (|24D and after some algebra, the first equality in (i) follows. To 
prove the second equality in (i), we differentiate 

(25) 2t(t -l)^ = -D' + D''' + tD^^ + ^ 

with respect to t along the characteristic t^/^e^''^^*^ = A, multiply the resulting equation 
by 2t{t — 1), and then substitute the t-derivatives by the corresponding right hand sides in 
Proposition ^A\ , the end result being 



(26) At{t - 1)'^ + 2(t - l)(3t -1)^-AD' = 1. 

Note that ( ^61) is the same as (^21) . Since the initial conditions D^{0,t) = Dl.{0,t) = are 
also the same as those for Z)°'^^, we conclude that the second equality in (i) holds. 
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We have, from Proposition |4.1| (iii) and (iv), that 

t(t _ + tD^' = {t- 1)^^ + D',^', 

CLL (JiL 

which can be written simply as 

(27) t|((t-l)Z^^^) = |((t-l)D°-^^). 

The right hand side of (pTf ) is readily computed, i.e.. 



dV^ ' ' dt \2{te^^(^-t) - 1) J dt\ 2(^2-1) 

-t-^/^A + 1 
2(^2-1) ' 

so that (127|) now reads 



Straightforward integration yields that 

2 2ri/2^ + lnt ^ f{A) 



D^'= — — + 



2(l-t)(l-A2) ' t-l 

2e-(i-*)+lnt ^ /(tV2e-(i-t)) 



2(l-t)(l-te2^(i-*)) t-l 
The condition D^'^{{),t) = then forces 

1 + lnt 



m 



l-t2 



With this /, and after some algebra, (ii) follows. 
We have, from Proposition |4.1| (i) and (ii), that 



(28) t{t - + D^'' =t\t-^)^+ tD\ 

By (i), 

_ t(e"(i-*) - 1)2 _ - 2t'/^A + t 
~ 2(t-l)(te2^(i-*) -1) ~ 2(t-l)(A2-l)' 
The right hand side of ( PBD is easily computed, i.e., 

2. JD^ r.^ At^/^-AH 



dt 2(t-l)(A2_i) 
so that (p8|) can be rewritten as 

+ 



dt t{t-l) 2(t-l)(A2-l)' 
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By the standard procedure for solving first order linear ODEs, we have that 
(29) 



t (2At'/^ - A^lnt 

t-i\ 2 - 1 ' 



2(l-t)(l-te2-(i-*)) t-1 
where g is an arbitrary function. The condition D^'^{0,t) = then implies that 

Reporting g back in p9|), and after some algebra, (iii) follows. □ 



Some properties of D\{t), etc., are readily deduced from Theorem |4l2 . 
Corollary 4.3. Forn ^ 2, 

(ii) DY{t)=t-D^\t-^), 

(iii) D^\t)=t''D'^\t-^), 

(iv) D'^>\t)=t-D'^^\t-'), 

(v) Dl{t)=rDl{t~^), 

Proof. Since D^{x, t) = D^'^'^{x, t), equating the coefficients of of D^{x, t) and D'^'^'^{x, t) 
(i) follows. Since 



D^'^ixt t'^\ - 

^ ' ^ (i_t-i)(i_t~ie2-t(i-t-^)) 

[_1 + _ t) + e-"(i-*)] 
(l-t)(l-te2^(i-*)) 

(l_t)(l_teMi-i)) ^ 

equating the coefficients of x^ on both sides, (ii) follows. To prove (iii), we replace t by 
and X by xt in D^'^{xt,t~^) = D^'^{x,t), the end result being 

D^''{x,t) = D^\{xt){t-^), (t-')-') = D'^'\xt,t-^), 

from which (iii) follows. By a calculation similar to that in (ii), we have that D'^'^'^{xt, t~^) = 
D^'^^{x,t) from which (iv) follows, (v) follows from (i) and (iv). □ 

Remark 4.4. Corollary [4.3| (ii)-(v) state that -D^(t), etc., have certain reflectional symme- 
try. For example, the polynomial t"'D^'^{t~^) is obtained by reflecting the coefficients of 
D^'^{t) about Thus, Corollary |J (ii)-(in) (resp., (iv)-(v)) say that D^'^{t) (resp., 

Dl^{t)) are obtained from D^'^{t) (resp., D^'^'^{t)) by reflection, and vice versa. Corollary [4.3| 
(i) further says that Dl^[t) and D^^'^{t) are themselves reflectionally symmetric. However, 
D^^{t) and D^'^it) do not enjoy this property. 
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Corollary [4.3| (i) enables a furthur simplification of the recurrence relation for Dnif) derived 
in the previous section. 



Corollary 4.5. We have 

Dn{t)= [(2n-l)t+lp„_i(t) + 2t(l-t)D;_,(t) 

+ {l-t)[t-'D'„,\{t)-tDf_,{t)]- 

Proof. By Corollary ^]3| (i) and (|13D, the last three terms on the right hand side of Corol- 
lary |3.4| can be written as 

- t)D'„,Ut) + (1 - t)'Df_,{t) + 2(1 - t)D'^^ht) 
= (1 - t)[{l + t-')D',:\it) + (1 - t)Df_,{t) + Di_,{t) + D'^^lit)] 
= (1 - t)a^.,{t) + (1 - t)r'D^_,{t) - tDr_,{t)]. 

Aftering regrouping terms, the corollary follows. □ 

Remark 4.6. The first line of is precisely the recurrence relation for the Eulerian 
polynomial -B„(t) of type B. 

Denote by D{x,t) the exponential generating function for the Eulerian polynomial Dn(t) 
of type D, i.e., 

(31) D{x,t) = Y,Dn{t)^, 
where Doit) = 1. Since -Di(t) = 1, it is clear that 

(32) D{x, t) = l + x + D\x,t) + D^^\x, t) + D^'^{x, t) + D°'^^(a;, t). 
Computing D{x,t) by ( |5^ and Theorem |4.2| , we obtain the following result of Brenti |]T|. 
Corollary 4.7. We have 

(1 - t)e"(i-*) - xt{l - t)e2-(^-*) 



(33) D{x,t) 



1 - te2^(i-*) 



We have now the exponential generating functions for the sub-Eulerian polynomials at our 
disposal. By reversing the procedure described in the paragraph following the definitions of 
D^{x,t), etc., we have the following results. 

Corollary 4.8. Forn ^ 2, 

^^^_ = ^|2"-i[(A; + i)« + fc"]-(2A;+l)"}t^- ^ 



t(l-t)"-i > J ^ i(l_i)n-r 



ii) 7^^^ = E[(2^ + 1)" - + 1)(2^)1^'. 



(m 



t) 



" \^_^ = V[(2/t + 1)" - T{k + 1)" + nT'-\k + l)"-i]t^ 
1 — tr 



tH\-t) 
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Proof. We only prove the first equality in (i); the remaining assertions follow from similar 
reasoning. 



n>2 



t(i-t)«-in! t h-r ' 

e - 1)2 



2(1 -te2^) 



2 



e 



= ^[2"-i(A; + 1)" - {2k + 1)" + 2"-iA;"] — . 

Equating the coefficient of x" on both sides, the result follows. □ 



5. A RECURRENCE RELATION 

In the previous section, we computed the PDEs which the generating functions of the sub- 
Eulerian polynomials satisfy. By solving these PDEs, we obtained the generating functions. 
It is important to note that the PDEs encode the recurrence relations for the sub-Eulerian 
polynomials. This procedure for computing generating functions can be reversed, i.e., start- 
ing from the generating function, we determine a PDE of minimal order which the generating 
function satisfies; the desired recurrence relation then follows upon equating coefficients. 

We shall apply this reversed procedure to -D(x, t) to obtain a recurrence relation for -D„(t). 
Proposition and ( p2D suggest the partial differential operator Hit — V)dt + (1 — 2xt)dx 
for D{x,t). The proofs of the following lemma are easy, albeit tedious, exercises of calculus 
which we omit. 

Lemma 5.1. We have 

(i) m - m + (1 - 2xt)dx]D{x, t) = ^ l_te2^a-t) — 

(ii) [2t{t - m + (1 - 2xt)d.]'D{x, t) = i^^^' + J ■ 



1 - te2^(i-*) 

(iii) [2t{t - l)dt + (1 - 2xt)dx]^ = At\l - tfdtt - 4t(l - t)(l - 2xt)dxt + (1 - 2xtfdxx 
2t(l - 2x)dx + 4t(l - t)(l - 2t)dt. 



Proposition 5.2. The function D = D{x,t) satisfies the following PDE: 

4t2(l - ty[x{l + t)- l]Du - 4t(l - t)(l - 2xt)[x{l + t) - l]Dxt+ 
(1 - 2xtf[x{l + t) - 1]D,, + 2t(l - t)2[-2 + (3 + t)x]Dt+ 
[At - x{l + Sty + 2xH{3 + 2t + 3t'^)]Dx + (1 - tfD = 0. 



16 CHAK-ON CHOW 

Proof. Use Lemma to form the expression 

a[2t{t - l)dt + (1 - 2xt)d^fD{x, t) + (3[2t{t - l)dt + (1 - 2xt)d:,]D{x, t) 
(34) _ (1 - t){[a{l + t)+ I3{1 + 3t^)]e"(i-*) + {-at - 2/?t^)e^"(^-*)} 

1 - te2^'(i-*) 

where a = a{x,t), and (3 = l3{x,t) are functions to be determined. Now assume that the 
right side of (^) is equal to 

(35, 

where 7 = 'j{x,t) is a function to be determined. Equating the coefficients of e^'-^"*^ and 
g2x(i-t) numerators of fBl) and f|35!), we have that 



a{l + t) + + =j, and - at - 2pt^ = -jxt, 

whose formal solutions are a = 7(1 — + 3t^) — 2t], and (3 = —7(1 — + t) — 1]. 

Setting 7 = (1 — t)^, we get a = [x(l + 3t^) — 2t] and /? = — + t) — 1]. Reporting these 
choices of a, /?, and 7 in (|3^), and using Lemma |0|(iii), the result follows. □ 



It is trivial to see that 

_ _ ^ ^ _ ^ / / \ ^ ^ T / \ 



We are now ready to state the main result of this section. 

Theorem 5.3. Forn ^ 1, the Eulerian polynomials Dn{t) of type D satisfy 

Dn+2{t) = [n{l + 5t) + At]Dr,+i{t) 

+At{l-t)D'^^,{t) 

+ [(1 - ty - n{l + 3tf - 4n{n - l)t(l + 2t)]Dn{t) 
-[4nt(l - t)(l + 3t) + 4t(l - t^D'^it) 
-4t\l-t)^D';{t) 

+ [2n{n - l)t(3 + 2t + 3t^) - An{n - l)(n - 2)t'^{l + t)]Dn-i{t) 
+ [2nt(l - t)2(3 + t) + 8n(n - 1)^2(1 - t)(l + t)]D'^_^{t) 
+Ant\l-tni + t)D';^,{t). 
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Proof. Denote by /, //, . . . , VI the successive terms of the PDE as in the previous proposi- 
tion. Then 



X 



X'' 



I = [-At'il - tf + At\l - tfil + t)x] J2 Kit) 

//= [4t(l-t)-4t(l-t)(l + 3t)x + 8^2(1 -t)(l + t)x2]^D;+i(t) 

/// = [-1 + (1 + 5t)x - 4t{l + 2t)x^ + 4t\l + J2 Dn+2{t) — , 

IV = [-4t(l - tf + 2t{l - tns + t)x] J2 D'nit):^, 

V = [4t - (1 + 3tfx + 2t(3 + 2t + 3t^)x^] ^ A^+l(t) 

VI = {l-trJ2Dnit)^, 



n>0 



n>0 



SO that 



nl [n — 1)! 

, At{l-t)D'^^,{t) 4t(l-t)(l + 3t)D;(t) 



n\ {n — 1)! 

U\l~t){l+t)D'^_^{t) D^+^it) , (l + 5t)D„+i(t) 



+ 



(n-2)! n\ (n - 1)! 

4t(l + 2t)Dn{t) ^ U\l + t)Dn-l{t) 



(n-2)! (n-3)! 
4t(l-t)2D;(t)^2t(l-t)2(3 + t)Z};_i(t) 



{n — 1)! 
4tD„+i(t) (1 + 3t)2Z}„(t) 



n\ {n — 1)! 

2t(3 + 2t + 3t2)D„_i(t) (l-t)2D„(t) 



+ 



(n-2)! n! 

where [x"] (/ + ■■■ + V^J) denotes the coefficient of a;"" in / + ■ ■ ■ + VI. Multiplying through 
by n\ and collecting terms, the theorem follows. □ 



The recurrence relation (^) involves Dn{t) and its derivatives only. By equating coeffi- 
cients of t'^, we have the following recurrence relation for Dn,k- 
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Corollary 5.4. 

Dn+2,k = {n + 4:k)Dn+l,k 

+ {5n-Ak + 8)Dn+i,k~i 

+ [(1 - n) - 4(1 + n)- Ak{k - 

+ [-2(1 + n + 2^2) + 8(1 -n){k-l) + 8{k - 1)(A; - 2)]Dn,k-i 
+ [(1 -n- 8^2) - 4(1 - 3n){k - 2) - 4(A; - 2){k - 3)]Dn,k~2 
+ [6nk + 4:nk{k - l)]Dn-i,k 

+ [6n(n - 1) + 2n(4n - 9){k - 1) - 4n(A; - l){k - 2)]D„_i,fe_i 

+ [An{n - 1)2 + 2n(fc - 2) - An{k - 2){k - 3)]D„_i,fc_2 

+ [2n(l - 3n + 2^2) + 2n(5 - 4n)(A; - 3) + 4n(A; - 3)(A; - 4)]Dn-i,k-3- 

Theorem f).3\ (or equivalently Corollary ^.4|) constitutes the solution to a problem con- 
cerning the recurrence relations for the Eulerian polynomials Dn{t) (or the Eulerian number 
Dn,k), posed by Brenti Note that some of the coefficients of Dn,k are not linear in n as 
well as in k. A direct combinatorial proof of Corollary |5.4| is desired, however. 

A closely related problem, also posed by Brenti, is whether the Eulerian polynomials Dn{t) 
have real zeros only. That 6„,(t) and -B„,(t) having only real zeros follow easily from their 
recurrence relations. Given the complicated nature of the recurrence relation for -D„(t), the 
corresponding assertion for -D„(t) need not follow in a similar manner. 
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